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Abstract 

In addition to the dielectric and magnetic permeability constants, another constant is generally 
needed to describe the electrodynamic properties of a linear isotropic medium. We discuss 
CN ' why the need for the third constant arises and what sort of physical situations can give rise to 

' a non-zero value for it. This additional constant, which we call the "Activity Constant" and 

■ denote by can explain optical activity and other phenomena from a purely macroscopic and 

. phenomenological point of view. 

^ ' 
(N 

a; 
a.: 

In principle, the best way to discuss electromagnetic processes taking place within a material medium 
[ is to use the subatomic picture in which matter is but a huge number of fundamental particles crowded 
in the vacuum. The charges and currents of these particles can be fed into the equations of vacuum 
^ I electrodynamics. The resulting equations, if solved, can describe all electromagnetic phenomena. 
■ - ' In practice, however, this is not a viable approach because, as already commented, the number of 

charges in any macroscopic sample of a material is enormous. One, therefore, has to take a simplified 
approach. This is achieved by dividing all the sources of charge into two parts: the external and the 
induced. The external charges are described by constructing the charge- and current- densities of them 
and feeding them into Maxwell's equations. The presence of the induced charges is taken into account by 
introducing a few empirical constants in the equations. 

The number of such constants necessary depends on the nature of the medium in question. For the 
simplest isotropic materials, it is well known that one needs two quantities, the dielectric constant e and 
the magnetic permeability fi, to describe all phenomena of macroscopic electrodynamics [^. However, 
for certain isotropic systems one needs in general three constants (not two) to describe electrodynamic 
phenomena in them. While this fact is also well known 0, it is not widely discussed in standard physics 
textbooks on electromagnetism Q. 

In this paper we explain, in simple terms, why the third constant (which we denote by C) arises and what 
kind of situations might lead to the vanishing of this extra constant. Some interesting physics associated 
with is discussed. This includes the phenomena of optical activity of some active solutions, and of atomic 

1 



Table 1: Electromagnetic quantities and their dimensions, in Heaviside-Lorentz (HL) units, in terms of the fundamental 
units of length {L), time (T) and electric charge (Q) To find the dimensions of the Fourier transforms of these 
quantities, multiply any entry by L'^T. In the last two columns, we give the recipe to transform the quantities in other 
systems of units. 



Quantity 


Dimension 
in HL units 


Replace by the following to change to 


Gaussian units 


SI units 


P 






p/V^ 


j 


QL-'^T-^ 


\/4^j 




E 




E/^/ivf 




B 


QL-^ 






c 




c 


l/V/^oeo 



parity violation which arises due to the breaking of the parity symmetry in the weak interactions. Because 
is crucial in explaining these phenomena from a macroscopic point of view, we refer to it as the "Activity 
constant". A discussion of the same issues from a quantum field theoretic point of view can be found in 
Ref. I] and Ref. |]. 

2 Formulation of Macroscopic Electrodynamics 

In the vacuum, Maxwell's equations consist of two homogenous equations in the electric field E and the 
magnetic field B: 

(9B 

V-B = , cVxE = -— (2.1) 

dt ^ ' 

and two inhomogeneous equations involving the sources: 

5E 

V-E = p , cVxB- — =j. (2.2) 

We have used the Heaviside-Lorentz units here so that no factors of Am appear in the equations. On the 
other hand, Coulomb law and Biot-Savart law appear now with factors of (47r)~^. As usual, c is a universal 
constant denoting the magnitude of the velocity of light in the vacuum. For rewriting any equation of the 
present paper in other popular systems of units, see Table |^. 

For future purposes, let us note the dimensions of various electromagnetic quantities in these units. Of 
course, p has the dimension of QL~^ since, by definition, it is the charge per unit volume Once that is 
known, Eqs. ( |2.1| ) and ( |2.2| ) determine the dimensions of E, B and j. The results are summarized in Table 

0- 

It is convenient to re-write all the above equations in terms of their Fourier components. We therefore 
express every quantity <I>(x, i) as a superposition of plane waves through the relation Q 

$(x, t) = j j ^^1^ ^(k' ^) exp[i(k • x - ujt)] . (2.3) 

From the usual results of Fourier analysis, it is then easy to show that 

V ^ ik, d/dt -iuj , (2.4) 
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which indicates how the derivatives on co-ordinate space functions affect the Fourier components. With 
the aid of these properties, we can cast Maxwell equations in the form 



k- B = 
ik-E = p 



ck X E = wB 
ick X B + iLoF, = j 



(2.5) 
(2.6) 



This form makes the ensuing discussion much easier since the equations become algebraic. This is especially 
convenient later when we introduce the dielectric constant e and the magnetic permeability of a medium 
through Eqs. ( 2.1C| ) and ( ^.13| ), which would become integral relations in co-ordinate space. Notice that 
we write the Fourier components with the same notation as the space-time functions themselves. This will 
not cause any confusion because the space-time functions will never again be used in this article. 

As mentioned in the Introduction, the discussion of macroscopic electrodynamics starts from the ob- 
servation that p and j contain, in addition to any external sources (denoted by the subscript "ext"), the 
sources induced in the medium (denoted by the subscript "ind" ) : 



Pcxt + Pind 



J — Jcxt ~l" Jind • 



(2.7) 



Ordinarily, the induced sources are expressed in terms of the polarization (P) and the magnetization (M) 
as follows: 



Pind = • P 

jind = —iLoF + ick X M . 



(2.8) 



In this form, the expressions for pi^^ and jind mimic those for total p and j given in Eq. ( |2.6D , and they 
automatically satisfy the equation of continuity, which takes the following form in Fourier space: 



k • jx - ^/Ox = , 



(2.9) 



where the subscript 'x' stands for the induced sources. This signifies the conservation of induced charges. 



The conservation equation (2^) also holds when 'x' stands for the total or just the external charges. 

The dielectric constant e and the magnetic permeability p for a system are now introduced by writing 



(6-l)E 



M 



[l-p-')B. 



(2.10) 



In writing these relations it is assumed, as we do in this paper, that the response of the system to applied 
electromagnetic fields is linear. This excludes special systems such as ferroelectrics and ferromagnetics 
from consideration, but otherwise does not impose any real restriction provided that the applied fields are 
not extremely large. Further, we will consider the simple case of isotropic media, for which e and p can be 
taken as scalars |8|. One can then use Eqs. (2^), (12.81) and (|2.10D to rewrite Eqs. (|2.5|) and (^.61) as 

k-B = 

ik • D = pext , 



^ and (|2.10D to rewrite Eqs. (|2^) and 

ck X E = a;B 

ick X H + ituY) = jext • 



where 



D = E + P = eE 



H = B - M 



-1 



B. 



(2.11) 
(2.12) 

(2.13) 



Since the source terms are just the external charges and currents, one can in principle solve Eq. (2.12) 
for D and H. Once the empirical constants e and p are known, one can then use Eq. ( p^ ) to derive E 
and B. 

In the next section, we will see that this treatment makes an important assumption tacitly, and because 
of that, it cannot describe an important class of physical phenomena within the framework of macroscopic 
electrodynamics . 
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3 The most general parametrization of medium properties 



3.1 Parametrization of the induced sources 

The crucial question in macroscopic electrodynamics is: how to parametrize the induced sources pind 
and jind- It is important to recognize that these quantities cannot be determined from considerations of 
macroscopic physics alone. In standard textbooks of electromagnetism, some simple microscopic models of 
the medium are invoked to express P and M in terms of microscopic quantities, and thereby motivate the 
parametrization discussed above. However, such models may not describe all possible media and, therefore, 
the expressions in Eq. ( |2.8| ) need not be the most general ones. From a purely macroscopic point of view, 
we can try to parametrize pind and jind in the most general form subject to the conditions of isotropy and 



linearity. As we show below, this approach indicates that Eq. ( p. 81) indeed does not provide the most 
general expressions for the induced charge and current densities. 



To clarify the point, we summarize Eqs. ( |2.8D and ( 2.1C| ) by writing 



And = -i(e-l)k-E (3.1) 
jind = -i(e-l)a;E + ic(l-p"^)kx B. (3.2) 

Looked at this way, it seems that we are trying to write down pind and jind as linear functions of the electric 
and the magnetic fields, where the coefficients depend on to and k. One can then ask, are the combinations 
the most general possible ones subject to the assumption of isotropy? 

The answer is 'y^s' for the expression for pind in Eq. (^]^). Note that p is a scalar as far as rotations 
are concerned, and the only scalars which are linear functions of E and B are k • E and k • B. Of them, 
the latter is identically zero (Eq. (|2.5|)), so that And has to be proportional to k • E, as shown in Eq. ( p.l] ). 

The situation is different for Eq. (13.21) . Here, jind is a vector, and there are four vectors involving E 
and B, viz., wE, ck x B, wB and ck x E, which have the same dimensions as j. The last two are equal 
by Eq. ( |2.5D , so that only the first three can be taken as independent. One can therefore write the most 
general expression for j as 

ji^d = -i(e-l)a;E + ic(l-p"^)kxB-iCwB. (3.3) 

Notice that the coefficient in the first term is fixed to be the same as the coefficient in Eq. (|3.1| ) by current 
conservation. Thus, in general, we have three dimensionless constants describing the electromagnetic 
properties of an isotropic medium: e, p and C,. 

Looking back at Eq. ( |3.2D , it now becomes immediately obvious why it is incomplete. Since jind is 
a vector having three components, one needs, in general, three basis vectors to write it down. In Eq. 
( |3.2| ), we used only two basis vectors, wE and k x B. Therefore, Eq. ( |3.2D is valid only under some extra 
assumptions, the precise nature of which will be discussed in the next section. 

One comment is in order. The last equation in Eq. (^]^) expresses the total current j as a combination 
of the two vectors wE and k x B. If we have had no complaint about that, why do we argue about Eq. 
I), which expresses jind as a linear combination of the same two vectors? 



The answer to this question lies in the difference of the nature of Eqs. ( p.6| ) and (^.2|). Eq. (2^) is a 
law, describing the fundamental nature of electromagnetic interactions. Eq. ( |3.2| ) is a parametrization of 
the induced current, which expresses our ignorance about the details of the dynamics of all processes going 
on in the medium or our inability to quantify them in minutest details. Such parametrization must be as 
general as possible subject to the conditions of the problem, and Eq. ( |3.3D provides this general form for 

Jind- 

With this parametrization of the induced sources, the pair of inhomogeneous equations given in Eq. 
(|2.6|) become 

iek • E = pext , 
4 



ieajB + i—k x B + {(ajB = jcxt , (3.4) 



while the sourceless equations ( p. 5]) remain intact. These two equations plus the equations in ( |2.5D can be 
taken as the fundamental equations for the E and B fields in a medium in the presence of external sources. 
In the next sections we will elucidate some of the physical implications and consequences associated with 
the constant 

3.2 Constitutive relations 

As already mentioned in the Introduction, the inadequacy of the simple constitutive relations given in Eq. 
( 2.13| ) for describing a certain class of physical phenomena has been recognized for a long time In many 



instances this deficiency has been addressed g |, |l^, 0] by postulating Eqs. (|]l|) and (|]l|) together 
with a more general set of constitutive relations to replace those in Eq. ( 2.13| ): 



D = eE + /3B, (3.5) 
H = 7E + ;U~^B. (3.6) 

With this approach, it seems that one needs two extra constants for the medium, /? and 7, rather than 
the constant ( discussed earlier. However, it is easy to see that the two constants introduced here are not 
independent. For this, substitute Eqs. ( |3.5| ) and ( ^I^ ) into Eq. ( p.l2| ). This yields 

iek • E = pext , 

i-k X B + ieujB + i{(3 + j)ujB = jext , (3.7) 

Thus, it is only the sum /3 + 7 that appears in the field equations, and therefore only three constants have 
an independent physical significance. Moreover, comparing these equations with those in Eq. (^.4|), we see 
that /3 + 7 = C. Because of this redundancy, one can try to put an arbitrary constraint to fix (3 and 7, e.g., 
one might put /3 = 7. 

While the approach of postulating more general constitutive relations rather than the simple ones given 
m Eq. (|2l|) has been followed in the literature for a long time, the lesson we have learned from the above 



discussion is that, in general, the constitutive relations that have been used contain redundant parameters. 
The reason for the redundancy is easy to understand. The componets of E and B are not independent. 
They are constrained by Eq. ( |2.5| ). Similarly, the components of D and H also are not independent. The 
constitutive relations, being relations between these two sets of variables none of which are independent. 



can also have redundant parameters. This can also be seen from Eq. ( 2.12 ). In fact, only the equation 
for jext is independent there; the equation for pext can be derived from it using the charge conservation 
equation, Eq. (|2.9| ). Taking the general expression for jcxt from Eq. (3^) and eliminating B using Faraday's 
law, Eq. ( |2.11|) , the equation for jext in Eq. (|2.12| ) becomes 

ck X H + wD = ewE H k x (k x E) + cCk x E (3.8) 

Thus, we are trying to define two vectors H and D from this single equation, in terms of just E, which is the 
only independent variable here. There is no unique way to do this, and any parametrization will contain 
some arbitrariness. In contrast, the components of jind are all independent. By focusing the attention on 
them, we have succeeded in finding a parametrization of the induced sources that is both minimal and 
complete. Because of this redundancy of the parametrization in terms of the constitutive relations, it is 
more convenient to use the parametrization of jind instead, which we will do in the rest of this article. 

We should comment that the restriction of the number of electromagnetic constants to three applies 
only in the Fourier space, which we use in the present article. In co-ordinate space on the other hand, 
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Table 2: The transformation properties of various physical quantities under the operations P, T and C. The quantity rjp, 
for example, is defined in Eq. ( |4.4[) , and the definitions of rjT, "qc etc are similar. The last two columns also show the 
transformation under the combined operations CP and CPT. 





Transformation given by 


— : 

Quantity 


VP 


Vt 


Vc 


Vcp 


VCPT 






+ 




+ 


+ 


B 


+ 
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k 
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+ 






u 


+ 




+ 


+ 




P 


+ 


+ 








j 








+ 




A 








+ 





one may add many terms using higher and higher derivatives of the E and B fields in the definition of the 



induced current [13|. Translated to the Fourier space, those additional terms are taken into account by 



allowing the constants e, n and C to depend on higher powers of k and lo. 

4 The constant ( and the discrete symmetries of space-time 

Why is it that the constant is not introduced for the simplest isotropic materials? Alternatively, what 
are the conditions under which ( can have a nonzero value? The key to the answer to these questions lies 
in the discrete symmetries of space and time, as we discuss in this section. 

4.1 Definition of the discrete symmetries P, T, C 

For orientation, consider Newton's law for a particle of mass m moving under the infiuence of the gravita- 
tional force exerted by another particle of mass M placed at the origin. The equation of motion is 

GMm 

ma = 2 — r, (4.1j 

where the right hand side denotes the gravitational force on the particle. If we perform space inversion 
(i.e., change r to — r), the force changes sign, but so does a, since it is defined to be cPr/dt"^. So Eq. ( [4.1D 
does not change. Similarly, if we perform a time inversion (i.e., change t to —t), neither side of Eq. ( [4. ID 
changes so, again, the equation remains invariant. 

For a long time, it was believed that these two discrete operations, viz. of space inversion or parity (P) 
and time reversal (T), do not change any laws of physics. If this were true, all the terms in a physically 
acceptable equation would transform the same way under P and T, as in the example discussed above. 



Indeed, we can verify that such is the case with the Maxwell equations in (2J) and (^^) or, equivalently, 
with their Fourier space counterparts given in Eqs. ( |2.5D and (|2.6| ). To show this, we proceed as follows. 

Consider, for example. Parity. The identification given in Eq. (^]^), implies that in Fourier space the 
parity operation is induced by the substitution 

(k,u;) ^(-k,u;) (4.2) 

It is very easy to verify that, if we simultanously make the following substitutions 

E(k,w) -E(k,u;), 
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B(k,cj) 


p ^ 


B(k,w), 




j(k,u;) 


p ^ 


-j(k,u;) , 




p{k,uj) 


p ^ 




(4.3) 



then the equations are not changed. If for every quantity, including oj and k themselves, we write the 
above transformation rules in the form 

$ (k, LJ) 7?p$ (k, LO) , (4.4) 

then the multiplicative constant r/p has the value given in Table |2|for the various quantities. 

The same procedure is followed for the Time reversal symmetry, for which Eq. (^^) suggests |jl^ 

(k,a;)^(k,-a;). (4.5) 

The fact that the Maxwell equations are invariant under P and T further bolstered the idea that all the 
laws of physics are invariant under those operations. 



It is easily seen that the Maxwell equations (Eqs. ( ^.l] ) and (2^), and their Fourier space version given 
in Eqs. (|2.5| ) and (p.6|)) have another discrete symmetry since they are unaffected by the operation under 
which p, j, E and B all change sign. This operation, called charge conjugation (C) because the sign of all 
the charges are reversed, plays an important role in quantum field theory, where it converts a particle into 
its associated antiparticle. 

4.2 Implications of P, T, C symmetries on e, yU, ( 

Now that we know the transformation properties of the electromagnetic fields and currents, we can look 
back at Eq. (|3.3| ) with respect to the discrete symmetries C, P and T, and some important combinations 
like CP and CPT. To proceed we first notice that we can consider instead 

jext(k, w) = ie(k,uj)uj'E,(k,uj) H -ik x B(k, cj) + iC(A:, u;)a;B(k, w), (4.6) 

p{k,uj) 



which is equivalent to Eq. (3.3) because of the definition in Eq. (2.7). We have explicitly indicated the 
dependence on k and oj of all the quantities involved, because it is important to keep that in mind in 
what follows. In addition, we have explicitly used the fact that the quantities e, p and C, which represent 
properties of the medium, can depend on k only through its magnitude = |k| since we are considering 
an isotropic medium. 

Before we proceed any further, we want to note that there is one relation that the electromagnetic 
constants e, /x and C must satisfy based on very general grounds. It is simply the requirement that the 
fields E, B and the current density j are, in coordinate space, real quantities. This implies that, for real 
values of uj, the electromagnetic constants satisfy 

r(-k,-u;) = /(k,u;), (4.7) 

where / stands for any of them. As noted earlier, the dependence on k can be only through its magnitude 
k since we are restricting the discussion to isotropic materials. Thus, in our case, we can write 

rik,-u;) = f{k,u;), (4.8) 

which means that the real parts of these electromagnetic constants are even in uj whereas the imaginary 
parts are odd. We stress that this is a very general property, not tied up in any way with the discrete 
symmetry operations that we are going to discuss now. The discrete symmetries, if they are valid in any 
given situation, impose extra conditions on the electromagnetic constants. 
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Consider, for example, Parity. The question we want to answer is the following: If we take each quantity 
that appears in Table ||, including k and uj, and make in Eq. ( |4.6| ) the replacement indicated in Eq. ( [4.4[ ), 
what relations must be satisified by e, (j, and C, in order for that equation to remain the same? In this 
fashion we easily derive the consequences of the parity symmetry on e, fj, and and similarly with the 
other discrete symmetries. We discuss them one by one. 



Parity : The parity reversed form of Eq. (4.6) is 



or 



- jg^t = -iewE - -ik X B + iQuo'B (4.9) 



jext = ^ewE + -ik X B - iC.uj'B . (4.10) 



Here and below, the arguments of e, /i, and unless otherwise given, should be assumed to be k and w. 
Comparing Eqs. ( [4. 61) and ( [4.10| ), we conclude that 

Parity symmetry (" = . (4-11) 



Charge Conjugation : The charge conjugated form of Eq. ( [4. 6]) is the same as the original one since 
each terms changes sign, and therefore we do not get any constraints on e, /x, and Q. 

CP : Obviously, this gives the same constraint as parity alone does, since charge conjugation does not 
affect Eq. (|4.6| ), as stated above. 



Time reversal : The time-reversed form of Eq. ([4. 61) is 



- jext = -iew (/c, -uj) E -ik X B + iCw [k, -uj) B, (4.12) 

^^{k,-u) 

which gives, on comparison with Eq. ([4.6|), the following constraints: 

e (/c, — w) = e{k,uj) 
Time reversal symmetry =^ ^ jjL{k,—uj) = fj,{k,uj) . (4-13) 

C(/c,-w) = -C{k,uj) 



Taken together with Eq. (|4.q ), this implies that ( must be purely imaginary and an odd function of oj, 
whereas e and fi should be real and even functions of uj. 

CPT : In this case, we get similarly the following constraints 

{e{k,—u)) = e{k,u;) 
H{k,-uj) = n{k,uj) . (4.14) 
ak,-uj) = Q{k,u) 

4.3 Constraints on e, /i, C in ordinary media 

From the relations deduced above, it is quite clear that the presence of the (^-term implies some proper- 
ties that are asymmetric under P and CP transformations. Similarly, since T and CPT imply mutually 
inconsistent constraints on C,, at least one of them must be violated in order for Q to exist. 
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We said earlier that classically it was believed that all interactions conserve P, T and C. From this, one 
might have expected C to vanish. However, we know now that P, T and C are all violated in fundamental 
interactions between particles. 

Since 1957, we know that parity is violated by weak interactions. Since 1964, we know that CP is also 
violated by the same, although to a much lesser extent. Since jmd arises because of complicated processes 
taking place within the medium, including weak interaction processes, there is no reason why the ^-term 
should not be present in Eq. 

If we leave the discussion at that point, it might seem that the parameter can at best be very small, 
since the weak interaction is indeed very feeble. What could be worse, the (^-term violates CP, and the 
CP-violating effects are known to be smaller than the usual weak interaction strength by about three orders 
of magnitude. So (, if present, would seem to be extremely tiny. But that need not be true in general, as 
we argue shortly. 

If we had carried the previous line of argument further, we would have concluded that the only constraint 
on the value of C is C (^i ^) = C (^) which follows from CPT symmetry, since there are very strong 
reasons to believe that fundamental interactions between particles can never violate CPT. This again, 
would have been wrong for the reason discussed below. 

In principle, the quantities e, fi, ( could be determined if we could solve the equations of motion of the 
particles that compose the medium. These equations include the effects of the mutual interactions between 
the particles as well as the effects of the local E and B fields. These are not external prescribed fields but 
are the fields that enter into Maxwell's equations, which must also be solved for, and which in turn contain 
the induced charge and current densities in the source terms. Thus, Maxwell's equations plus the equations 
of motion of the particles in the medium form a closed system which must be solved self-consistently. To 
date, no one has succeeded in solving this formidable problem and in practice one turns to a more modest 
approach. 

In classical as well as quantum physics, the calculation of e, /i, C proceeds in two steps. In the first, 
the equation of motion of a single particle is set up, which includes the effects of the local fields, as well 
as the effects of the other particles in the medium. Prom this, the contribution of a single particle to the 
induced charge and current is obtained. In the second step, this single-particle contribution is suitably 
averaged over the particles in the medium, thereby obtaining an expression for the induced charge and 
current densities, from which the values for e, fx, C can be extracted. 

It is now clear that the asymmetries with respect to any operation can enter into any quantity in either 
of these two steps. In the first step, it enters if the underlying interactions do not respect the symmetry 
generated by the operation, as is the case of weak interactions with respect to the symmetries of P, C and 
CP. But the asymmetries can also creep into the process of taking the average over the particles in the 
medium, because the medium might have a surplus of one kind of particles over others. 

As a well-known example, consider the mass of the electron and its antiparticle, the positron. In 
vacuum, the two masses must be equal. Formally, this is ensured by the combined symmetry CPT. Inside 
a chunk of normal matter, however, there are many many electrons and hardly any positrons, so that 
the interactions of the two particles become very different. As a result, they move with different effective 
masses within a medium. Thus, the CPT breaking effects creep in through the asymmetry in the medium. 

In the case of our present interest, considerations of CPT and CP symmetries are irrelevant for the 
same reason. In fact, the CP and CPT asymmetries are large since normal matter does not contain any 
antiparticles at all, and hence neither CP nor CPT stand in the way of obtaining a large C. 

Let us now remind ourselves that in order to obtain a large we also need a large P-asymmetry. The 
P-symmetry is of course violated by weak interactions, and that can give rise to a small value of To 
obtain a large we need a large P-asymmetry in the medium. 

Take a free electron gas, for example |15|. An electron in this gas can be specified by its momentum p 
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and spin angular momentum s. Under the action of parity, p — p, s ^ s. However, in the metal, there 
are as many electrons going in the +x direction, say, as in the —x direction . So the parity operation does 
not affect the overall momentum distribution of the metallic electrons. In other words, this medium does 
not have any P-asymmetry and so C must arise in this medium through parity violation in weak interaction. 

Similarly, under the T operation, p — p, s — s. But again, since the distributions of both p and 
s were random to start with, the overall distributions do not change by this operation, i.e., the medium 
is T-symmetric. Therefore, barring the small T-violating effects which are 10~^ times weaker than the 
P-violating effects in weak interactions, the consequences of T-symmetry should be taken seriously. As we 
found before in Eq. ( [4.13| ), this implies that C must be purely imaginary and an odd function of uj. 

This is quite striking. Note that e and fi are both even functions of oj, so that it is conceivable to think 
of a medium where they are constants. In this case the medium is non-dispersive. But C cannot have a 
nonzero constant term and therefore all media are dispersive when it comes to C- Of course, we emphasize 
that this comment is true only to the extent that T-violation in fundamental interactions can be neglected. 

When we consider more complicated media, it is possible to obtain P-asymmetry from the medium itself. 
A well-known example is that of helical molecules, with one helicity in excess over the other. This occurs 
in natural sugar solutions, for example, where in fact one helicity is absent altogether. If we consider the 
passage of electromagnetic waves through such a solution in which the helical axes are randomly oriented, 
and if the wavelength is larger than the molecular dimension so that the wave cannot "see" individual 
molecules, we have an isotropic medium with a built-in P-asymmetry. This gives rise to large ^-term, and 
physically it manifests itself as the phenomenon of optical activity, as we discuss in Sec. ^ 

In short, there is no reason why ^ should be zero in normal matter, and its value can be large if 
there exist large asymmetries in the medium. Therefore it is imperative that we understand its physical 
significance. 



5 Dispersion relations 

It is well-known that a great deal about the nature of wave-propagation is learned from the dispersion 
relations. These are the relations that co and k must satisfy in order that solutions to Eq. ( |3.4D exist when 
/9ext = and jext = 0; i.e., 

ek-E = (5.1) 
eu;E + ■^kxB + CujB = 0. (5.2) 

These two, plus the two equations in Eq. (p.5|), can be taken as the fundamental equations governing the 
E and B fields in a medium in the absence of external charges and currents. 

In the vacuum, since e = ^ = 1 and C = Oi Eq. ( ^ ) clearly shows that E must be perpendicular to 
the direction of propagation of the electromagnetic wave. Thus, one can discuss electromagnetic waves in 
terms of two basis polarization vectors, ei and 62, both being perpendicular to k, as shown in Fig. |l]. 

In the medium, Eqs. ( |5.lD and ( p-2\ j can be satisfied for a longitudinal electric field if B = and u 
is such that e = 0. Since the physical properties of this longitudinal mode depend on the expression for 
e only, we focus our attention on the transverse modes in order to discuss the physical implications of 
Thus, eliminating B from Eq. ( |5.2D by the use of Faraday's law and using k • E = 0, we obtain 

uj'^B = —k^E - —Lokk X B , (5.3) 
efi e 

where k = k/A;. Let us first see what happens if = 0, which is the case treated in standard textbooks 



10 



k 




Figure |I| 

Figure 1: The relative orientations of the two transverse polarization vectors with respect to the wave vector. 



[16|. In this case, Eq. (|5.3D immediately gives the dispersion relation 



en 



(5.4) 



In non-dispersive media (i.e. when e and fj, are independent of uj and k), this gives the speed of propagation 
as cj yfe[i. Apart from the transversality condition k • E = 0, there is no restriction on E. Therefore, 
we conclude that the two transverse modes have the same physical characteristics and, in fact, any linear 
combination of them is equally appropriate for discussing the transverse wave properties. 

In general, however, the two transverse modes have different properties when ^ does not vanish. To 
proceed in this case, we take the cross product of Eq. (|5.3| ) with k and use the transversality condition 
• E = to obtain 



uJ^k X E 



(? - cC 
—k^k X E + — w/cE. 
e/x e 



Taking linear combinations, this pair of coupled equations for E and A; x E can be rewritten as 



Therefore, solutions exist only if lo and E satisfy 



e/i eke. 



Ezbi/c X E 



e/i eke / 
E = iXkxB, 



(5.5) 

(5.6) 

(5.7) 
(5.8) 



with A = ±1. Thus, the two transverse modes no longer have the same dispersion relation. Furthermore, an 
arbitrary linear combination of them does not, in general, obey the wave equation (|5.3| ). The combinations 
that travel as waves are the ones that satisfy Eq. ( ^.8]) , which are in fact the positive and negative circular 
polarization vectors e ± ife x e, where e is any vector transverse to A: [p 



We argued before, in connection with Eq. (4.13), that e and fi are real whereas C is purely imaginary 
if time-reversal invariance is satisfied. In this case, it follows from Eq. (5.7) that that cu is real for real 
values of k if \(^\ < ke\fi\ /uj. The amplitude of a monochromatic wave will thus remain unchanged during 
the propagation, which means that there is no absorption in the medium. However, if 7^ 0, the speeds of 
propagation of the two circlularly polarized transverse modes are different. 

In a more general case, time-reversal symmetry can be violated. Then, for a given k, uj is complex in 
general, so that the exponential factors in Eq. ( |2.3D give some damping (or growth) with time. Physically, 
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Figure |^ 

Figure 2: The effect of a Parity transformation on the wave vector k and the transverse polarization vectors. The empty 
circles in diagrams (a) and (c) indicate that the direction of ei is towards the reader, whereas the filled circle in diagram 
(b) indicates the direction to be into the paper. 



this corresponds to absorption by (or emission from) the medium. Even in the case of C = 0, absorption 
of the two transverse modes can arise if e or ^ is complex. But the absorption has to be equal for the two 
modes. For C / 0; that is not the case. 

For the sake of simplicity, we will stick to an unabsorbing medium. This means that u) is real for all 



wave vectors k. From Eq. (5/7), this implies that e, /i, are real, C, is purely imaginary and \C,\ < kc\^\ /cj. 
As already remarked, in this simple case the effect of a non-zero value of C is to produce a difference in 
speed between the two circularly polarized degrees of freedom. 



6 Physical implications of ( 

We now discuss the conclusions of the last section in a more intuitive way on the basis of the discrete 
operations P, C and T. Again, we concentrate on the transverse modes only, and draw the polarization 
vectors ei, 62 along with the unit wave- vector k in Fig. ^a. 

Now imagine the effect of parity on the vectors in Fig. Surely k reverses direction, since it is a polar 
vector denoting the direction of wave propagation. To see the effect of parity on ei and 62, we first need a 
precise definition of them. This depends on the way we choose the scalar and vector potentials to express 
the electric and magnetic fields. For the sake of convenience, we take the "radiation gauge" in which the 
scalar potential vanishes and 



E = iujA, B = -ik X A . 
In this gauge, we can express the vector potential as 

A = ^ aAekA , 

A 



(6.1) 



(6.2) 



where the quantities a\ are numerical coefficients in classical electromagnetism. Now, under parity, E 
reverses sign, so that from Eq. ( |6.1D , we see that A must do the same. From Eq. ( |6.2D , we then obtain 
that the vectors must also change sign. All of these facts have been summarized in Fig. 
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However, since we are dealing with isotropic systems, the absolute directions of the vectors do not 
matter, only their relative orientations do. Thus, rotating Fig. ^ around the e2-axis by 180°, we obtain 
Fig. Comparing this with Fig. we see that effectively, for the same k, the effect of the parity 
transformation has been to reverse the sign of 62- 

An alternative way to discuss the polarization of electromagnetic waves is to use, instead of the linear 
polarization vectors ei and 62, the right and left circular polarization vectors defined by 

e± = -^{ei ± 162) ■ (6.3) 



Since 62 = k x ei, these are of the form ii ±ik x ei, as mentioned after Eq. (5.8). In this language, thus 



the effect of parity is easily understood. Since effectively 62 — > —62, as shown earlier, we conclude that 
under parity 

e+ < — > e_ . (6.4) 

Thus, if parity is a good symmetry, there should be no difference in the physical properties of a right- 
circular and a left-circular polarization state of an electromagnetic wave. Such is the case, for example, in 
the vacuum, since Eqs. (|2.5| ) and ( p. 61) are invariant under parity. 

In matter, as we pointed out, this is not the case. Eq. ( |3.3D and its corollary in Eq. ( |4.6| ) break parity 
because of the (^-term. Thus, in a medium, the right- and the left- circularly polarized electromagnetic 
waves can have different physical properties. For example, they travel with different speeds through matter, 
have different indices of refraction and coefficients of absorption. 

It is easy to see that this difference manifests itself as optical activity. Consider a plane polarized wave 
incident on a medium at x = 0, moving in the positive x-direction. Suppose the direction of polarization 
of the incident wave is along ip = 0, where ip denotes the azimuthal angle in the plane perpendicular 
to the direction of propagation. As is well-known |l^, this plane polarized wave can be treated as the 
superposition of a right-circularly polarized wave with phase angle <f+{0) = u)t and a left-circular polarized 
wave of the same amplitude with phase angle ip- (0) = —ijot. The zero in the paretheses here denote the 
fact that these phases are given for the interface where x = 0. 

Now the wave propagates inside the medium, where C 7^ 0, and consequently, the same value of u) 
corresponds to different values of the wave number k for the two circularly polarized components, as we 
argued before. Let us call these values of k by /c+ and k^. Then, after traversing a distance x within the 
medium, the phase angle for the right- and left- circularly polarized components become 

ipj^{x) = ujt — kj^x, ip-{x) = —ujt + k-X , (6.5) 

so that the resultant direction is given by 

1 1 

fix) = - [f+{x) + ip-{x)] = - (A;_ - k+) X . (6.6) 

This shows that the direction of polarization of the wave changes with x within the medium, which is 
precisely the phenomenon of optical activity. 

A similar treatment can be applied to the phenomenon that occurs in the optical-rotation experiments 
that have been performed to test the parity-violating weak force between the electrons and nucleons. 



which coexists with the electromagnetic force [19|. In these experiments, a beam of lineraly polarized light 
traverses a cell containing atomic vapor of bismuth, and the direction of polarization of the wave changes as 
the beam traverses the cell. From a macroscopic point of view, the effect can be understood and described 
as discussed above. 

The last two equations, together with our earlier analysis leading to Eq. ( |5.7| ), bring out clearly the 
fact that the difference between A;_ and k+ is introduced because C / Iii view of this, the quantity 
( may be called the activity constant of a medium. 
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7 Conclusions and summary 



We have explained why, in order to describe electromagnetic phenomena in an isotropic medium, one 
needs in general three parameters, or electromagnetic constants. Two of them are introduced and studied 
in detail in every textbook on Electromagnetic Theory, and are called the dielectric constant and the 
magnetic permeability of the medium in question. The third one, which we denote by the symbol (" in this 
paper, is not as widely discussed or not introduced at all in physics textbooks. We showed that ( has some 
features which are very different from the other two. For example, a medium can have a non-zero value of 
( only if there are asymmetries of parity, CP etc. These asymmetries may come either from a violation of 
those symmetries in fundamental interactions between particles, or from the constituents of the medium. 
Since it is known that particle interactions violate parity and CP, the value of C is in principle non-zero for 
all normal media. Morever, if the constituents of the medium are parity-asymmetric objects like helical 
molecules, the effect can be large. 

Physically, the presence of the C term introduces differences between the characteristics of the right- and 
left-circular ly polarized electromagnetic waves. These differences manifest themselves in various physical 
phenomena, e.g. optical activity and atomic parity violation experiments. Our approach in the present 
paper provides a general framework to study and interpret these phenomena from a unified point of view. 



Appendix 



A Time reversal of the Maxwell equations in Classical and Quantum 
mechanics 

As we have discussed in the text, a non-zero value of C has deep implications regarding the discrete space- 
time symmetries of the system. The transformation rule that we have used for Time- Reversal is apropriate 
in the context of Classical Physics which, for clarity, we have adopted in the text. In the realm of 
Quantum Mechanics that definition of Time-Reversal cannot be used. Therefore, it is important to discuss 
the appropriate transformation rule in Quantum Mechanics, since in principle C has to be calculated in 
that context. In this appendix, we discuss the different definitions of the time-reversal transformation and 
why they arise. 

Let us use as starting point Eq. ( |4.6| ) which gives the external current density in terms of the electric 
and magnetic fields. For the sake of convenience, we define the following quantities: 

Lf(k,u;) = -ie'^'^^k^ + iLoC5''^ , (A.l) 

where the superscripts a, h denote Cartesian components, 6°"^ denotes the Kronecker delta symbol and e""^^ 
is the completely antisymmetric tensor with e^'^^ = +1. With these, we can rewrite Eq. (4^) as 



jcxt = E^l''E'' + LtB^ (A.2) 

b 

Since jext and B change sign under time reversal whereas E does not, we can say that the Time- Reversal 

invariance of the Maxwell equations is equivalent to the statement that Le and Lb have the following 
properties: 

L£;(k, -w) = -L£;(k,Cj) 

LB(k,-a;) = LB(k,u;). (A.3) 
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These in turn imply the relations on e, and C that have been obtained in the text. This definition of 
the Time-reversal operation can be called the classical definition. If Eq. (A. 3) holds, then the Maxwell 
equations are invariant under the substitutions 



E(k, uj) E'(k', io') = E(k, io) , 

B(k, Lo) B'(k', uj') = -B(k, uj) , 

j(k,u;)^ j'(k',.;') ^-j(k,..), (A.4) 

where 

(k,a;)^ {k',Lo') =(k,-w). (A.5) 



Eq. ( A.4 ) is, of course, equivalent to the identification of the value of t^t for the various quantities given 
in Table 

In Quantum Mechanics, however, these transformation rules cannot be used. The reason is that, while 
they leave the Maxwell equations invariant, these rules do not leave invariant the quantum mechanical 
equations of motion for the charges, namely the Schrodinger equation. A simple way to understand this is 
to recall that, in Quantum Mechanics, one identifies u = £ /fi and k = p/^, where £ and p denote the energy 
and momentum of the particle and fi is Planck's constant divided by 27r. Now, from the non-relativistic 
definitions of momentum and kinetic energy of a point particle [viz., p = and £\^i^ = |m ], one 
expects that p changes sign under time-reversal whereas £ does not. Thus, instead of Eq. ( [A.5| ), the 
transformation is defined by 

(k,u;)^ (k',cj') =(-k,a;). (A.6) 

This is how the transformation must be defined in order to leave the Schrodinger equation for a free 
particle invariant. The interesting point is that this transformation is a symmetry operation of the Maxwell 
equations also, provided that Le and Lb satisfy, instead of Eq. ( |A.3D , the following: 

L^(-k,u;) = -L£;(k,a;), 

L|j(-k,a;) = Lb{Koj). (A.7) 

If these relations hold, then the Maxwell equation is invariant under the substitutions: 

E(k,tc;)^ E'(k',cj') =E*(k,a;) 

B(k,cu)^ B'(k',u;') =-B*(k,u;) 

j(k,u;)^ j'(k',^') =-j*(k,u;) 

i ^ i' =i* = -i. (A.8) 

We will refer to these transformation rules as the Quantum Mechanical ones. Since the transformation 
rules are different from the ones discused before, the consequences of time reversal symmetry on e, /x and 
change as well. 

The last comment is meant to emphasize that Eqs. (A.4) and ([A.8|) are two different transformation 



rules, and the symmetries that they imply are also different. However, they are equivalent if, but only if, 
Le and do not have an absorptive part. In that case, those quantities depend on i only through the 
combinations ik and itu, and therefore taking the complex conjugate of them is the same as changing the 
sign of u> and k in their argument. It is then immediately seen that Eqs. ( |A.3| ) and ( |A.7| ) are really the 
same. However, if L^; and Lb have an absorptive part, then the two transformation laws are different. 
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The classical transformation is a symmetry operation only if Eq. ( A.3| ) holds. That equation implies 
the relations which have been derived in Eq. ( 4.13| ), which in turn imply that e, fi are real while C, is purely 
imaginary; i.e., neither have an absorptive part. The classical transformation law is not a symmetry oper- 
ation if any of those quantities has an absorptive part. However, the quantum mechanical transformation 
law is a symmetry operation even when they have an absorptive part; what is needed is that Eq. (A. 7) 
be satisfied, and that is always true because e, /i and C, satisfy the relation ( |4.g| ). As already discused, in 
classical physics this relation is a consequence of the fact that the fields and the current are real quantities, 
while in Quantum Mechanics it is just the requirement of Hermiticity of the current and the fields. 

We summarize the situation as follows. In a non-absorbing medium, the classical and quantum defini- 
tions of the Time-Reversal operation coincide. However, the classical Time-Reversal symmetry is broken 
if e, /i or have an absorptive part. Therefore, in the realm of Classical Physics, it is impossible to get 
an absorptive part of those quantities except if the theory breaks the (classical) Time-Reversal symme- 
try. The Schrodinger equation breaks that symmetry; quantum effects break that transformtion. But the 
equations of Quantum Mechanics (and of course Maxwell equations) are invariant under a more general 
transformation, which continues to hold even when e, and C have an absorptive part. The existence of 
the C, term implies that P is broken but not T (with the quantum mechanical definition of T) nor C. Thus, 
C breaks P, CP and CPT. 
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